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Abstract 

Let A be a real symmetric matrix of size N such that the number 
of the non-zero entries in each row is polylogarithmic in N and the 
positions and the values of these entries are specified by an efficiently 
computable function. We consider the problem of estimating an arbi- 
trary diagonal entry {A"^)jj of the matrix A™ up to an error of eb"^, 
where 6 is an a priori given upper bound on the norm of A, m and e 
are polylogarithmic and inverse polylogarithmic in TV, respectively. 

We show that this problem is BQP-complete. It can be solved effi- 
ciently on a quantum computer by repeatedly applying measurements 
of A to the jth basis vector and raising the outcome to the mth power. 
Conversely, every quantum circuit that solves a problem in BQP can 
be encoded into a sparse matrix such that some basis vector |j) cor- 
responding to the input induces two different spectral measures de- 
pending on whether the input is accepted or not. These measures can 
be distinguished by estimating the mth statistical moment for some 
appropriately chosen m, i.e., by the jth diagonal entry of A^. The 
problem is still in BQP when generalized to off-diagonal entries and it 
remains BQP-hard if A has only —1, 0, and 1 as entries. 

*e-mail; janzing@ira.uka.de 
^e-mail: wocjan@cs.ucf.edu 



1 Introduction 



It is believed that a quantum computer is more powerful than a classical 
computer in the sense that it makes possible to obtain faster algorithms 
for some computational problems than the best classical algorithms. How- 
ever, it is still not understood well enough which problems are tractable for 
quantum computers. It is therefore be desirable to better understand the 
class of problems which can be solved efficiently on a quantum computer. 
In quantum complexity theory, this class is referred to as BQP. Meanwhile, 
some characterizations of BQP are known ^ [21 El 15- Here we present a 
new characterization of BQP which is related to the computation of powers 
of large matrices. 

It should not be too surprising that computational problems can be for- 
mulated in terms of "large" matrices. For example, the transformations of a 
quantum computer can be represented by multiplication of matrices of a cer- 
tain type. However, the matrix problems derived from this representation, 
would usually not be very natural in classical terms. They are, of course, 
natural, as physical questions about the behavior of quantum systems. For 
instance, the problem of estimating the entries of products of unitary ma- 
trices which are given by a tensor embedding of low-dimensional unitaries, 
is BQP-complete, but it is not obvious where problems of this nature could 
arise in real-life applications referring to the macroscopic world. 

It is known that Hamiltonians with finite range interactions can generate 
sufficiently complex dynamics that can serve as autonomous programmable 
quantum computers jH]. Therefore, it is not unexpected that problems re- 
lated to spectra and eigenvectors of self-adjoint operators lead to computa- 
tionally hard problems. One could think that many of such problems could 
be solved efficiently on a quantum computer. However, results proving that 
questions pertaining to the minimal eigenvalue of Hamiltonians are QMA- 
complete (HI El EI demonstrate that it is unlikely to find efficient algorithms 
for this problem. 

The situation changes dramatically when we do not aim at deciding 
whether some Hamiltonian H has an eigenvalue below a certain bound 
but only whether a given state \ip) has a considerable component in the 
eigenspace corresponding to a particular eigenvalue of H. This problem can 
be answered by (1) applying i?-measurements to \tp) several times and (2) 
statistically evaluating the obtained samples. It has been shown in Q that 
measurements of so-called /c-local operators^, applied to a basis state, solve 

^An operator on n qubits is called fe-local if it can be decomposed into a sum of terms 



2 



all problems in BQP. This proves that some class of problems concerning 
the spectral measure of fc-local self-adjoint operators associated with a given 
state, characterize the class of problems that can be solved efficiently on a 
quantum computer. Unfortunately, the requirement of /c-locality restricts 
the applicability of these results since it is not clear where fc-local matrices 
occur apart from in the study of quantum systems. For this reason we have 
constructed a problem with sparse matrices that does not require such a 
fc-local structure and show that a very natural problem, namely the compu- 
tation of diagonal entries of their powers, characterizes the complexity class 
BQP. 

The paper is organized as follows. In Section |2l we define formally the 
problem of estimating diagonal entries. In Section we prove that this 
problem can be solved efficiently on a quantum computer. To this end, we 
use the quantum phase estimation algorithm to implement a measurement 
of the observable defined by the sparse matrix. To do this it is necessary 
that the time evolution defined by the sparse matrix can be implement effi- 
ciently. Since the diagonal entries of the mth powers are the mth statistical 
moments of the spectral measure, we can estimate them after polynomi- 
ally many measurements provided that the accuracy is sufficiently high. In 
Section |3 we show that diagonal entry estimation encompasses BQP. The 
proof relies on an encoding of the quantum circuit which solves the consid- 
ered computational problem into a sparse self-adjoint matrix such that the 
spectral measure (and hence an appropriately chosen statistical moment) 
corresponding to the initial state depends on the solution. In Section El 
we generalize this result to matrices with entries —1, 0, and 1. The idea 
is that the gates, which are encoded into the constructed Hamiltonian are 
not required to be unitary, even though the circuit that then realizes the 
corresponding measurement is certainly unitary. This fact could be inter- 
esting in its own right. For example, it could be possible that there are even 
more general ways of simulating non-unitary circuits by encoding them into 
self-adjoint operators. In this context, it would be interesting to clarify the 
relation to other measurement based schemes of computation [HI IIUI ITT] . 

2 Definition of diagonal entry estimation 

Before we define the decision problem "diagonal entry estimation" we in- 
troduce the notion of sparse matrices and spectral measures. Here we call 
an N X N matrix A sparse if it has no more than s = polylog(A^) non-zero 

that act on at most k qubits 0) 
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entries in each row and there is an efficiently computable function / that 
specifies for a given row the non-zero entries and their positions (compare 

[IlllSlIIil). 

Let ^4 be a self-adjoint matrix of size N x N and 

A 

be the spectral decomposition of A. Let \^) be some normalized vector 
of size A^. The spectral measure induced by A and is a probability 
distribution on the spectrum of A such that the eigenvalue A occurs with 
probability ||(5a|V')P- I^i the sequel we will repeatedly make use of the 
following observation. The expectation value of A*" in the state I?/') is given 

by 

{^\A"-m = Y,^"'mxm, 

A 

i.e., by the mth statistical moment of the spectral measure. 

In ^ , eigenvalue sampling is defined to be a quantum process that allows 
us to sample from a probability distribution that coincides with the spectral 
measure induced by A and \ip). Throughout the paper we refer to such a 
procedure as measuring the observable A in the state We now state the 
considered problem in a formal way. 

Definition 1 (Diagonal Entry Estimation) 

Given a sparse real symmetric matrix A of size N , an integer j € {!,..., N}, 
and a positive integer m = polylog(A/'), estimate the diagonal entry {A'^)jj 
in the following sense: 
Decide if either 

(A™),, ><7 + e6™ 

or 

{A^)j,<g-eh^, 

for given g G 6"*] and e = 1 /poly log (A^), where b is an a priori known 

upper bound on the operator norm of A. 

Problems of this kind arise, for example, in graph theory. Let A be the 
adjacency matrix of a graph with N vertices and degree bounded from above 
by s. Then the diagonal entry {A"^)jj of the mth power of A is equal to the 
number of paths of length m that start and end at the vertex j. 

It is important to note that the scale on which the estimation has rea- 
sonable precision is given by b"^. If the a priori known bound on the norm 
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is, for instance, 6' := 2b instead of b, then the accuracy is changed by the 
exponential factor 2™. Our results show that quantum computation out- 
performs classical computation in estimating the diagonal entries (provided 
that BQPt^ BPP). But one has to be very careful on which scale this result 
remains true. 

3 Diagonal entry estimation is in BQP 

To show that diagonal entry estimation is in BQP we briefly recall the formal 
definition of this complexity class |6|. 

Definition 2 (The class BQP) 

A language L is in BQP if and only if there is a uniformly generated family 
of quantum circuits acting on poly(r) qubits that decide if a string x of 
length r is contained in L in the following sense: 

y^|x,0) = ax,o|0) 8) |^x,o> +ax,i|l> ® l^x,i> (1) 

such that 

1- lox.iP > 2/3 z/x e L and 

2. |ax,iP < 1/3 i/x L. 

Equation ^ has to be read as follows. The input string x determines the 
first r bits. Furthermore, I additional ancilla bits are initialized to 0. After 
Yr has been applied we interpret the first qubit as the relevant output and 
the remaining r + I — 1 output values are irrelevant. The size of the ancilla 
register is polynomial in r. 

We now describe how to construct a circuit that solves diagonal entry 
estimation. Without loss of generality we may assume 6 = 1 and rescale 
the measurement results later. The main idea is as follows. We measure 
the observable A in the state \j) and raise the outcome value to the mth 
power. The average value of the obtained values over large sampling con- 
verges to the desired entry. The measurement is done by (a) considering 
j4 as a Hamiltonian of a quantum system and simulating the corresponding 
dynamics Ut = exp{—iAt) and (b) applying the phase estimation algorithm 
to Uf. The proof that this works follows from a careful analysis of possi- 
ble error sources. These are (1) errors due to the statistical nature of the 
phase estimation algorithm, (2) statistical errors due to estimation of the 
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expectation value from the empirical mean, and (3) errors caused due to the 
imperfect simulation of the Hamiltonian time evolution. We show that all 
these errors can be made sufficiently small with polynomial resources only. 

(1) We embed A into the Hilbert space of n qubits, where n = [log2 A^] • 
Let us first assume that the unitary U := exp(i^) can be implemented 
exactly. We apply the phase estimation procedure which works as follows 
[15j . We start by adding p ancillas to the qubits on which U acts. The idea 
is to control the implementation of the 2'th power of U by the Ith control 
qubit. More precisely, we have the controlled gates 

t^;:=|0)(0|«®l + |l)(l|(')®C/2\ 

where the superscript (/) indicates that the projectors |0)(0| and |1)(1| act 
on the Ith control qubit, respectively. Note that the decomposition of Wi 
into elementary gates is obtained by replacing each elementary gate in the 
circuit implementing U with a corresponding controlled gate. Similarly, Wi 
is realized by applying the quantum circuit implementing the corresponding 
controlled [/-gate 2' times. Set W := W1W2 ■ ■ ■ Wp. The phase estimation 
circuit consists of applying Hadamard gates on all control qubits, the circuit 
W, and the inverse Fourier transform on the control qubits. The desired 
value a is obtained by measuring the control qubits in the computational 
basis. Let {tp) be an arbitrary eigenvector of U with unknown eigenvalue 
g«27rv5 £qj, gQjj^g phase 93 S [0, 1). In order to achieve that the phase estimation 
algorithm outputs a random value a G {0, . . . , 2^ — 1} such that 

Pr(|v9-a/2f| <r/)>l-0, (2) 

for some 9, r] > it is sufficient JH] to set 

p:= riog(l/r?)l + [log (2 +(1/(20))]. 

Let I ijj) be an eigenvector of A with unknown eigenvalue A G [—1,1]. 
In order to determine A approximately using the outcome a in a phase 
estimation with U = exp(i^) we proceed as follows. First, we have to take 
into account that > 1/2 corresponds to negative values A. Second, we have 
to consider that the scaling differs by the factor 27r. Finally, we may use the 
additional information that not all A in [— vr, vr) are possible, but only those 
in [—1,1]. All outputs that would actually correspond to eigenvalues A in 
[l,7r] and [— vr,— 1) are therefore interpreted as ±1, repectively. Therefore, 
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we compute values z from the output a by 



z 



' a{2-K/2'P) for < a < 2P/(27r) 

1 for 2P/(27r) < a < 2^/2 

-1 for 2P/2<a<2P-2P/(2^) 

_ a(27r/2P) - 27r for 2^ - 2P/(27r) < a < 2^ 



This defines the random variable Z whose values z are approximations for 
A that satisfy the following error bound: 



This bound follows from ineq. by appropriate rescaling (note that our 
reinterpretation of values in [— tt,— 1] and [l,vr) explained above can only 
decrease the error unless it was already greater than vr — 1). Consequently, 
we have for every eigenstate of A with eigenvalue Aj the statement 



first term on the right hand side corresponds to the unlikely case that the 
measurement outcome deviates by more than 27r7/ from the true value. Since 
we do not have outcomes z smaller than —1 or greater than 1 the maximal 
error is at most 2. This leads to the error term 26. The second term 
corresponds to the case |Aj — z\ < 27rr], which implies |A™ — z™| < (27r7/)m 
because \i,z G [— 

We make the error in eq. © smaller than e /3 by choosing the parameters 
6 and r] such that 6 < e/12 and rj < e/(12 tt m). The number of control qubits 
can be chosen to be 



Pr(|A - Z\< 27rr/) > 1 - . 




p := 2[log((48m)/e)] . 



(4) 



This is sufficient since 



riog(l/r?)l + [log (2 + (1/(20)) < 2 [log ((48m)/e)l . 



(5) 



We decompose \j) into f/-eigenstates 



and obtain the statement 
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by linearity arguments and 

i 

Using the triangle inequality and the fact that the right hand side of ineq. Q 
is smaller than e/3 for each i we obtain 

\E\^^{Z^)-{A^),,\<e/^. (6) 

(2) Now we sample the measurement k times in order to estimate the 
expectation value Ey^{Z"^). Since we will later also consider the simulation 
error we want now estimate (A"^)jj up to an error of 2e/3. To achieve this, 
it is sufficient to estimate E\j^{Z^) up to an error of e/3. 

Let denote the average value of the random variable Z"^ after sam- 
pling k times. Since the values of Z^ are between —1 and 1 we can give an 
upper bound for the probability to obtain an average being not e/3-close to 
the expectation value. By Hoeffding's inequality ^1 Theorem 2] we get 

Pr(|Z^-i?l,)(Z-)| > I) < 2exp(^fc). 

In summary, we have shown for 6 = 1 that we can distinguish between 
the two cases in Definition ^ with exponentially small error probability. For 
6 7^ 1 we have to rescale the inaccuracy of the estimation by 6™. The whole 
procedure including repeated measurements and averaging can certainly be 
performed by a single quantum circuit in the sense of Definition [2 

(3) We now take into account that U = exp(zy4) cannot be implemented 
exactly. It is known that the dynamics generated by A can be simulated 
efficiently if A is sparse |121 1131 I14| . More precisely, for each t we can 
construct a circuit V which is 5-close to Ut = exp(— ^^4*) with respect to 
the operator norm such that the required number of gates increases only 
polynomially in the parameters n, t, and 1/5. We analyze the error resulting 
from using V instead of f/, where ||y — C/|| < 5. 

The phase estimation contains 2^"*"^ — 1 copies of the controlled-^ gate. 
Therefore the circuit Fy implementing the phase estimating procedure with 
V instead of U deviates from Fjj by at most 2^^^ 5 with respect to the 
operator norm, that is, \\Fu — Fv\\ < 2^+-"^ 6. 

Let q and q denote the probability distributions of outcomes when mea- 
suring the control register after the phase estimation procedure has been 
implemented with V and U, respectively. The li- distance between q and q 
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is then defined by 

Ik - q\\i ■= X] 1^*^") ~ 

ae{0,...,2P-l} 

where q{a), respective q{a), denote the probabilities of obtaining the out- 
come o. To upper bound \\q — q\\i we define a function s by s{a) := 1 if 
q{a) > a and s(a) := —1 otherwise. Let Q be the observable defined by 
measuring the ancillas and applying s to the outcome a. Then we can write 
||g — as a difference of expectation values: 

{^1;\fIqFu - FlQFv\tP) <2\\Fu - Fv\\\\Q\\ <2P-^^ 6 . 

This implies that the corresponding expectation values of Z"^ can differ by 
at most 5 because Z takes only values in the interval [—1, 1]. We choose 
the simulation accuracy such that 6 = e/{3- 2^^^) and obtain an additional 
error term of at most e/3 in ineq. Using that we have chosen p as in 
eq. (jH) we obtain that 6 G O(e^m^). 

Putting everything together we obtain a total error of at most e. Fur- 
thermore, this can be achieved by using time and space resources which are 
polynomial in n, m, and 1/e. This completes the proof that diagonal entry 
estimation is in BQP. 

It should be mentioned that off-diagonal entries {A'^)ij can also be es- 
timated efficiently on the same scale using superpositions \i) it \j) since the 
values {i\A'^\j) can be expressed in terms of differences of the statistical mo- 
ments of the spectral measure induced by those states. The scale on which 
the estimation can be done efficiently is then also given by e b"^ with an ap- 
propriately modified e which is still inverse polynomial in n. However, since 
BQP hardness requires only diagonal entries we have focused our attention 
on the latter. 

4 Diagonal entry estimation is BPQ-hard 

Now we assume that we are given a quantum circuit Yr that is able to decide 
whether a string x is in the given language L in the sense of Definition |2 
Using Yr we construct a self-adjoint operator A such that the corresponding 
spectral measure induced by an appropriate initial state depends on whether 
x is in L or not. Note that the proofs for QMA-completeness of eigenvalue 
problems for Hamiltonians have already used the idea to construct a self- 
adjoint operator whose spectral properties encode a given quantum circuit 
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Figure 1: Circuit U constructed from the original circuit Yr. Whenever the answer 
of the BQP problem is no, the output state of U is close to the input sate |x, 0) = 
\xi, . . . , Xn, 0, . . . , 0). Otherwise, the state |x, 0) is only restored after applying U 
twice. 

[HI 13 El- In these constructions, the existence of eigenvalues of a given 
Hamiltonian depends on whether or not an input state exists that is accepted 
by a certain circuit. In BQP, the problem is only to decide whether a given 
state is accepted and not whether such a state exists. Likewise, the problem 
is not to decide whether an eigenvalue of the constructed observable exists 
which lies in a certain interval. Instead, it refers only to the the spectral 
measure induced by a given state. This difference changes the complexity 
from QMA to BQP. 

The idea for our construction is therefore rather based on |1] which shows 
the BQP-hardness of approximate fc-local measurements. This result relied 
on the ideas in |E1 where the PSPACE-hardness of so-called exact fc-local 
measurements was proved. 

However, our description below will only at one point refer to these 
results since the observable we construct here is only required to be sparse, 
in contrast to the A;-locality assumed in |^ EI • In some analogy to 
we construct a circuit U that is obtained from Yr as follows: First apply the 
circuit Yr- Apply then a u^-gate. Implement then 1^^. The resulting circuit 
U is shown in Fig. ^ We denote the dimension of the Hilbert space U acts 
on by N. 

Let U be generated by a concatenation of the M elementary gates 
Uq, . . . , Um-1- We assume furthermore that M is odd, which is automati- 
cally satisfied if we decompose Yr in analogy to Yr and implement a cr^-gate 
between Yr and Yr . We define the unitary 

A/-1 

W:=Y, \l + l){l\(l?)Ui, (7) 

1=0 

acting on C^. Here the + sign in the index has always to be read 
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modulo M. We obtain 



A/-1 

=J2 \i){i\®Ui+M ■■■Ui+iUi. 

1=0 

Due to = 1 we have (W^^)'^ = 1. Thus, can only have the eigenval- 
ues ±1. This defines a decomposition of the space into symmetric 
and antisymmetric VF-invariant subspaces and S~, respectively with 
projections 

Q± :=i(l±VF^O- 
In the following we use the definition 

|sx) :=|0)®|x,0) 

for the initial state and restrict the attention to the span of the orbit 

jiy'lsx)} with Z EN. (8) 

Moreover, we use the abbreviations oq = ax,o and ai = ax,i- We consider 
first the two extreme cases =0 and |ai| = 1. If = the orbit Q 
is M-periodic and the action of W is isomorphic to the action of a cyclic 
shift in M dimensions, i.e., the mapping \l) i— > |(/ + 1) mod M), where \l) 
corresponds to VF'|sx) with / = 0, 1, . . . , M — 1. 

If |q!i I = 1 the action of W corresponds to a cyclic shift with an additional 
phase —1, i.e., the mapping |/) i— > |/ + 1) for ? = 0, 1, . . . , M — 2 and \M — 
1) 1-^ — 10). In the first case, the state |sx) induces a spectral measure 
being the uniform distribution on the Mth roots of unity, i.e., the values 
exp(— ZTT 2//M) for / = 0, . . . ,M — 1. In the second case, |sx) induces the 
measure R^^^ being the uniform distribution on the values exp(— ivr {21 + 
1)/M) for Z = 0, . . . ,M - 1. We observe that R^^^ and coincide up to 
a reflection of the real axis in the complex plane. 

In the general case, the orbit defines an 2M-dimensional space whose 
orthonormal basis vectors are obtained by renormalizing the vectors 

W^Q+\s^) and W^Q-\s^) with / = 0, 1, . . . , M - 1 . 

We obtain then a convex sum of R^^^ and R^^^ as spectral measures induced 
by W and |sx). The following calculation shows that |aop and |aip define 
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the corresponding weights: 

(Sx|g+|Sx) = ^(Sx|l + VF*^|5x> 

= i(x,0|l + C/|0,x) 

= i(i + (x,o|y>,y,|o,x)) 

= |aoP, 

where the last equahty follows easily by replacing l^.|x, 0) and its adjoint 
with the expression in eq. and its adjoint. Thus, we obtain the spectral 
measure 

We define the self-adjoint operator 

A := ^{W + W^). 

The support of the spectral measure corresponding to A is directly given by 
the real part of the support of R. To obtain the corresponding probabilities 
one has to take into account that in many cases two different eigenvalues of 
W lead to the same eigenvalue of A. 

To calculate the distribution of outcomes for ^-measurements we observe 
that R^^^ leads therefore to a distribution P^^^ on the (M — l)/2 eigenvalues 

Af)=cos^ for/ = 0,...,(M-l)/2 

with probabilities P^^'^ := 1/M and := 2/M for / > 1. Likewise, -R(^) 
leads to a distribution P^^^ on the (M — l)/2 values 

A«=cos^^^^ forZ = 0,...,(M-l)/2 

with probabilities i^(Aj_i)/2 = 1/^ and p/^^ = 2/M for / < {M - l)/2. As 
it was true for R^^^ and , the measures P^^^ and P^^^ coincide up to a 
reflection. 

We now set \j) := |sx), i.e., the input state is considered as the jth basis 
vector of C*^ C^. Then the diagonal entry {A"^)jj coincides with the mth 
moment of the spectral measure: 

(yl-),, = 0-|A-|i)=5^A-P(A), 
A 
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where A runs over all eigenvalues of the restriction of A to the smallest A- 
invariant subspace containing \ and P{X) denotes its probability according 
to the spectral measure corresponding to A. Since the latter is a convex sum 
of and P(^) we may write {A'^)jj as the convex sum 

(Ann = (l-|a,ftE(^r)"Vr + l«:I^E(\"')""^'" 

I I 

=: (1- |ai|2)£:o + |ai|2£;i. (9) 



The values £"0 and Ei can be considered as the mth statistical moments of 
random variables on [—1, 1] whose distributions are given by P(°) and P^^^ 
respectively. 

In order to see how the value {A^)jj changes with |q;i| we observe 

(M-l)/2 

£0- 2^ [\ ) ^ >Po +[\m-1)/2) - m + [\m-1)/2) ■ 
1=0 

Here we have used Aq'^'' = 1 and that the eigenvalues are numbered in a 
decreasing order. Thus, X(m-i)/2 is the smallest one. Due to the reflection 
symmetry of the measures we have Ei = —Eq. Now we choose m sufficiently 
large such that the term (^(^_i)/2)™ negligible compared to 1/M since 
we have then Eq — Ei k, 2 /M which is a sufficient difference for our purpose. 
In order to achieve this we set m := M^. We have 

2 4 2 

^S2-i)/2 = - cos(VM) >-l + ^- ^>-l + ^, 
where the last inequality holds for sufficiently large M. Due to 

^2 „ „2 



we conclude that 
and hence 



hm (1 - -^)^^ = 



(cos(7r/M))^' < (e-^)^ 



£^0>^-(e-^)^>7T7> (10) 



1 , 

(e 4 ^ 

where we have, again, assumed M to be sufficiently large. To see how {^A 
changes with |q;i| we recall 



jj 



{A^)jj = (1 - \a,\^)Eo + \ai\^E, = (1 - 2\ai\^) Eq , 



13 



by eq. © and the reflection symmetry. Using the worst cases jaip = 1/3 
for X L and |aip = 2/3 for x £ L we obtain 

and 

Using £^0 > 3/(4M) from ineq. (|T?H) we obtain 

1 



if the answer is no and 



otherwise. Setting g := (see Definition^ we may define e := 1/(4M). 
Then the diagonal entry is greater than g + e ii x L and smaller than g — e 
otherwise. The construction of A as the real part of a unitary ensures that 
ll^ll ^ 1 ='■ b. This shows that we can find an inverse polynomial accuracy 
e such that an estimation of the diagonal entry up to an error e b"^ allows to 
check whether x is in L. 



5 Generalization to matrices with entries 0, ±1 

So far we have allowed general real- valued matrix entries. We may strengthen 
the result of the preceding section in the sense that diagonal entry estima- 
tion remains BQP-hard if we allow matrix entries to be only 0, ±1. It is 
known that Toffoli gates and Hadamard gates form together a universal set 
for quantum computation We may thus replace the whole sequence 

Ui, . . . , Un of gates used in the definition of W (see eq. (|7I)) by a set of gates 
that consist only of TofFoli and Hadamard gates. Let T and H, denote the 
set of Toffoli gates and the set of Hadamard gates on C", respectively. We 
modify then the universal set and consider T^e/t UTright Ui?, where we have 
defined 7]e/t := TH and Tj-ight '■= HT. In words, TJe/t is, for instance, the 
set of gates that are obtained by applying an arbitrary Toffoli-gate followed 
by a Hadamard gate on an arbitrary qubit. The Toffoli gates are permu- 
tation matrices, whereas the Hadamard matrices have only entries ibl/\/2- 
Thus, the gates in T^jght and Tigjj have only entries ±l/-v/2, too. Therefore, 
all gates in Ti^ft U Tj-ight U H have only entries ibl/-\/2. Hence, the matrix 
A would only consist of such entries when using only gates that are taken 
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from our modified set of gates. By rescaling with \/2 we obtain a matrix A 
with entries 0, ±1. The rescaling is clearly irrelevant for the diagonal entry 
estimation problem since we have now spectral values within the interval 
[-V2,V2] and the accuracy required by Definition ^ changes by the factor 
(\/2)™ accordingly. 

6 Conclusions 

We have shown that the estimation of diagonal entries of powers of symmet- 
ric sparse matrices is BQP complete when the demanded accuracy scales 
appropriately with the powers of the operator norm. 

The quantum algorithm proposed here for solving this problem uses the 
fact that measurements of the corresponding observable allow to obtain 
enough information on the probability measure defined by the eigenvec- 
tor decomposition of the considered basis state. Given the assumption that 
BQP 7^ BPP, i.e., that a quantum computer is more powerful than a clas- 
sical computer, the required information on the spectral measure cannot 
be obtained by any efficient classical algorithm. This is remarkable since 
the determination of spectral measures is a problem whose relevance is not 
restricted to applications in quantum theory only. 
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